RADIO SCIENCE Journal of Research NBS/USNC-URSI 
Vol. 68D, No. 11, November 1964 



Geometrical Optics Convergence Coefficient for Whistler 

Propagation 

G. McK. Allcock 

Contribution From the Department of Scientific and Industrial Research, Physics and Engineering Laboratory, 

Lower Hutt, New Zealand 

(Received May 25, 1964) 

It is shown that the formula, given by Crary [1964] for the convergence coefficient for 
whistler signals propagated by multiple reflection between a convex earth and a concave 
ionosphere, can be considerably simplified. The simplified formula shows that the con- 
vergence coefficient for a given number of reflections is independent of the height of the 
ionosphere, and that antipodal focusing would be expected. The convergence takes place 
entirely in a horizontal direction transverse to the plane of propagation. 

A very simple first-order approximation to the formula is found to be adequate for most 
practical calculations. A second-order approximation shows that the convergence coeficient 
for a given propagation distance is almost ind< pendent of the number of reflections. 



1. Introduction 

Crary [1964] has recently derived an expression 
for the convergence coefficient in the ray theory 
equations for the field strength near the surface of 
a sharply bounded ionosphere due to a radio trans- 
mitter on the ground. His paper has served a 
useful purpose in showing that the convergence 
coefficient is indeed close to unity. This, as Crary 
points out, would intuitively be expected, since 
there are an equal number of reflections from the 
concave ionosphere and the convex ground. How- 
ever, equation (17) of (Vary's paper, which he has 
used for his illustrative calculations, is by no means 
in its simplest form, which we here derive. 

In figure 1 we have redrawn the first half-hop 
portion of Crary's figure 2, labeling the ground-based 
transmitter position G, the point of incidence of the 
ray on the ionosphere /, and the center of the earth 
0. We now construct OZ perpendicular to IG, 
ZY perpendicular to OG, and IX perpendicular 
to OG. Then 



OX=(a+h) cos (^) 



Also, since /OZY=/ZG v =<l>, 



OY=OZ sin 4>=a sin 2 <f>. 



(1) 



(2) 



Combining (1) and (2), 



(a+h) cos (2^7)-^ sin 2 <l>=OX-OY=XY 



=IZ cos <t> 

= (a+h) cos 6 n cos <£. (3) 
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Figure 1. Geometry for simplification of convergence coefficient 
formula. 



Also, 



a sin (j)=OZ=(a-\-h) sin B n . 



(4) 



We can now substitute (3) and (4) in Crary's 
equation (17), which reads 



otuh- 



(a+h) sin d n 



(2^+1) sin (j> tan </> sin ( • J (a+h) cos ( J— a sin 2 <£ 



(a+A) sin 7 cos 6 n 



1 J (27i+l) sin 2 sin ( y ) (a+h) cosfl w cos 4> I 

^ pn« A\( n.-\--h\ Gin 'v pn« # 



which reduces to 



^ = f (2n+l)sin( 2 -^) y / 
i. sin 7 J 



(5) 



This formula can also be derived analytically, though 
more laboriously, from Crary's equation (17), by 
applying the sine and cosine rules to the triangle 
IGO, and thence expressing the functions of (f> and 

d n in terms of functions of 



\2n+l/ 



2. Approximate Formulas for Distances Less 
Than 6000 km 

We can further simplify (5) for use at short dis- 
tances. Upon reference to figure 1, it can be seen 
that, for the case of VLF propagation where h is 90 

km or less, cos ( . j will always be greater than 
( ,7 ), except when propagation is tangential to the 
earth's surface. Therefore ( ) will always be 



\2n+] 

less than 0.168 radian (9.6°), 
approximation 



sin 



\2n+lf 



(2rc+l) 



will always be correct to within 0.5 percent, 
can also write 

7 3 
S1I17«7-- 



and thus the 



(6) 
We 

(7) 



to 1 percent accuracy for y less than 1 radian, i.e., 
for great-circle distances D less than 6370 km. 
Using approximations (6) and (7), (5) becomes 



1 



1/2 



1- 



= 1 + 



12 



'1 + 



WD\> 
12 W 



(8) 



The nature of the approximations is such that (8) 
will be correct to within about 1 percent for distances 
less than 6000 km. 



cos 4>{a-\-h) sin y cos 6 n 



It is seen from (8) that, to a first approximation, 
the convergence coefficient is independent of n, the 
number of reflections. A close estimate of the varia- 
tion of a wh with n is obtained by taking a second 
approximation 



sin 



2n- 



-1/ \2n+lJ G\2n+lJ 



(9) 



instead of (6), and proceeding as before, 
we o;et 



This time 



aa,»~l-+ 



n(n-\-l)y 2 
3(2ft+l) 2 



+ *»+;>. fgy 



3(2/i+l) 2 Va/ 



(10) 



which is equivalent to (8) for n— « . The error in 
using this approximation is less than 0.1 percent for 
distances up to 4000 km, and less than 0.5 percent 
at 6000 km. The dependence of a wh on n is ex- 
pressed by the term 0/rt --7- 9 ., which is almost in- 
F J 3(2?i+l) 2 

variant, having the value 0.07407 for n=l, and 
increasing slowly to 0.08333 for ti=°°. The sole 
exception is the case where n = 0, i.e., the transmitter 
and receiver are in direct line-of-sight. Here, as 
would be expected from elementary considerations, 

it— — — — =0, anda«A=l exactly. 

3(2n + l) 2 

The results of calculations of a wh for distances 
up to 6000 km are shown in figure 2. This shows 
graphically the virtual independence of a ah on n. 

3. Discussion 

One feature which is shown by (5) , but which was 
not evident in Crary's equation (17), is that the 
convergence coefficient in the whistler case is inde- 
pendent of the height of the ionosphere, except 
insofar as the minimum value of n is determined 
by the tangency condition <£=90. For a daytime 
ionospheric height of 70 km, the maximum valaes of 
D for n = 0, 1, and 2 are marked by vertical stops 
in figure 2. 

Equation (5) can also be derived ab initio by 
considering the geometry of a spreading bundle of 
rays with rectangular cross section, as used by 
Rawer [1952]. This approach shows that the net 
convergence in the whistler case takes place entirely 
in a horizontal direction transverse to the plane of 
propagation. By contrast, for the case of ground- 
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Figure 2. Convergence coefficient as a function of propagation distance. 



to-ground transmission, it is seen that a term similar 
to (5) appears in the formula for the convergence 
coefficient. [Crary's equation (1).] The remainder 
of the terms in the ground-to-ground formula, viz, 



\ a) L(l+h/a) cos (D/2an)-lJ 



state the amount of convergence taking place within 
the plane of propagation, in a direction perpendicular 
to the ray path. 

An examination of (5) shows that a uh becomes 
infinite when sin 7 = 0, provided that sin 

(- — -Wo, i.e., infinities occur for y = k-K (k = l, 

\2n+l/ 

2, 3 . . .). The principal value t = tt represents 
the case of antipodal focusing, already well known 
for ground-to-ground propagation. 



Finally, it should be emphasized that the above 
discussion has been confined to geometrical optics. 
No account has been taken of diffraction near grazing 
incidence, nor of other factors which enter into the 
calculation of field strength, such as the reflection 
coefficients of the ionosphere and the earth's surface, 
and the polar diagram of the transmitting antenna. 
This aspect of the problem has been discussed 
thoroughly by Wait [1961]. 
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